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Abstract. We describe results for the confinement-deconfinement phase transition as predicted 
by the Nambu-Jona-Lasinio model where the local and quantum Polyakov loop is coupled to 
the constituent quarks in a minimal way (PNJL). We observe that the leading correlation of two 
Polyakov loops describes the chiral transition accurately. The effects of the current quark mass on 
the transition are also analysed. 



The simultaneous ocurrence of two phase transitions where chiral symmetry is re- 
stored and hadrons become deconfined seems a unique and misterious feature of QCD 
matter at finite temperature (for an early review see e.g. |1]). Besides direct QCD lat- 
tice simulations, there exist theoretical constraints below and above the deconfinement 
phase transition. At low temperatures the leading thermal excitations correspond to a gas 
of weakly interacting pions |2D. Moreover, in the large N c limit with the temperature T 
kept fixed, if a chiral phase transition takes place it should be first order y|]. The use of 
resonance hadron Lagrangians implies that thermal corrections are 1 /N c suppressed 0]. 
At high temperatures one has a weakly interacting quark-gluon plasma (for a review 
see e.g. o]). However, the previous powerful constraints assume from the start a given 
phase and do not provide a clue on how chiral symmetry restoration and deconfinement 
are intertwined. 

The coupling of relevant order parameters such as the quark condensate for chiral 
symmetry breaking and the Polyakov loop for deconfinement at finite temperature can 
be made explicit in Polyakov chiral quark models, an amalgamate of colour and flavour 
degrees of freedom where the simultaneous chiral-deconfinement crossover can b e q uan- 
titatively studied with an acceptable phenomenological success |0, 0, HI El Ell EH EH 
El El EH]. Actually, in our recent work [10] we have shown why and how Polyakov 
loops must be coupled to Chiral Quark Models (CQM) to comply with large gauge in- 
variance at finite temperature and how the quantum and local nature of the Polyakov 
loop generates a rather sharp crossover at about the observed critical temperature, al- 
though uncertainties are expected. More specifically, ChPT and large N c constraints are 
naturally accomodated lll3ll within those models, hence solving a long standing puzzle 
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which was ignored for a long time; traditional CQM did produce a chiral phase transition 
while violating those restrictions. An immediate consequence of this new coupling is an 
upward shift of the critical temperature referred to as Polyakov cooling in Ref. llQfl . 

In this work we will deal with the PNJL model for definiteness. After bosonization 
the NJL Lagrangian reads 

r Q [S,P} = -Tr\og^+M + S + iy 5 P^ + ^rJ d 4 xtr f (S 2 + P 2 ) . (1) 

We use Tr for the full functional trace, tr^ for the trace in flavour space, and tr c for 
the trace in colour space. The UV divergencies in Eq. ^ from the Dirac determinant 
only affect in practice the zero temperature contributions II 1611 . In the Polyakov gauge 
Q = e lAA / T , where A4 is time independent and diagonal, and the minimal coupling is 
made d$ — > d$ — 1A4. Integrating further over the A4 gluon field in a gauge invariant 
manner ifivll yields a generic partition function of the form 

Z = J DSDPDQe- rG[il] e- r Q [S ' p > a] , (2) 

where DQ. is the Haar measure of the SU(A^ C ) colour group, Tq is the effective gluon 
action and Tq stands for the quark effective action. In general Q.(x) is a local and 
quantum variable since the gluon field itself depends on the point. As argued in ifioll 
mean field approximations 10 [ECU EH] generate a spurious gauge orbit dependence and 
a possibly complex Q. (violating colour charge conjugation) and they necessarily imply a 
non-vanishing value of the Polyakov loop in the adjoint representation, in contradiction 
with lattice results iflSll . 

For a constant value of the Polyakov loop, Q., and the scalar field S = M (which we 
identify with the constituent quark mass) the quark effective action is given by 



; T Q (M,a,T) = —Tr f (M-Mo) 2 -2N f j 



N c £ k 



V ^ ' ' ' AGs J " ' J J (2n)~ 

+TTr c jlog (l + e- Ek,T Q)j + log (l + e'^rf) } 



(3) 



where we have only retained the vacuum contribution, so there is no contribution of 
meson fields (S,P) (see II 1 311 for a chiral expansion up to 0(p 4 )), V is three dimensional 
volume and £^ = +Vk 2 +M 2 is the energy of a constituent quark with mass M. We 
define the Polyakov-loop averaged action 

e -r Q (M,T) = J dae -r G [n] e -r Q (M,a,T) (4) 



The value of M is determined by minimization of Tq(M, T) with respect to M, 
<9rg(M, r) /dM = 0, which corresponds to computing the integration in DSDP at the 
mean field level and determines M at a given temperature T, denoted as M* = M(T). In 



addition, the relation between the (single flavour) chiral quark condensate, (qq), and the 
constituent quark mass, reads 

2G s N f (qq)* = -(M*-m q ). (5) 

Any observable is obtained by using M* and averaging over £1. The integral in dQ. in 
the case N c = 3 and in the Polyakov gauge was computed numerically in Ref. [ 10]. Here 
we show a much simpler method which is based on evaluating the integral analytically 
for any N c in the low temperature limit and corresponds to take Q. small. (To see this use 
the formula / d 3 ke~ £k ^ T = AtzM 2 TK2 (y) , where Ki(x) is the modified Bessel function 
and K2(x) ~ ^Tt/2xe~ x for x —>■ °°). From Eq. © we get 

r Q (M,a,T) = r 2 (M,0)+2^ / £^^/^^e-'^/ r Tr c [a' 1 + ^. (6) 

n=l n ^ 

Expanding the exponent in Eq. © one obtains a power series in terms of Q. and The 
simplest correlation of two Polyakov loops is taken to be ifioll 

J dQTr c Q(x)Tr c Q?(y) = e - a ^ T , (7) 

with o the string tension, and yields the leading thermal contribution to the effective 
action 

T -T Q (M, T) = T -Y Q {M, 0) - TV a (lN f J ^- £k/T ) 2 + ■ " " , (8) 

where V a = SnT 3 /o 3 is the equivalent confinement correlation volume. As we see, 
the effect of quantum corrections on the Polyakov loop lowers the vacuum energy as it 
should. Moreover, they are 1 /N c suppressed, as one would expect in Chiral Perturbation 
Theory or a resonance gas model but unlike traditional chiral quark models without 
Polyakov loop. 

Minimizing with respect to the mass we get the effective temperature dependent 
mass M* and from Eq. © the corresponding (qq)* condensate can be evaluated. The 
approximate result is presented and compared to the full result [10] in Fig. [Hand, as 
we can see the approximation is quite efficient and very easy to implement in standard 
chiral quark models. For the Polyakov loop expectation value similar manipulations 
hold, yielding the leading order contribution 

/ 1 \ N f M 2 T . . , NfV a lM 3 T 5 m/t 

The full result and the approximated formula are compared in Fig. [l] In this case the 
agreement is only up to temperatures about 0.1 5To- 

The analysis above is done with physical current quark masses, m q = 5.5 MeV where 
one obtains T x = Tq = 256(1) MeV. This remarkable coincidence between transitions is 
not accidental nor depends on the particular choice of m q as can be seen in Fig.[U where 
we show the temperature dependence of (qq)* and L for m q — 0, 5.5, 40, 80, 120 and 
300 MeV. The corresponding susceptibilities are displayed in Fig. El 
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FIGURE 1. Chiral condensate (qq)* and Polyakov loop L = (ti c Q.) /N c . Left: Leading Polyakov loop 
correlation approximation (see Eqs. and @). We take the 2-flavor PNJL model, and ^fa = 425 MeV, 
f n = 93 MeV, M = 300 MeV, m u = m c / = m q = 5.5MeV. Right: Current quark mass dependence. 




FIGURE 2. Temperature dependence of dL/dT (left) and d(qq)*/dT (right) for several values of 
the current quark mass. The transition temperatures are T x = To = 255(1), 256(1), 266(1) MeV for 
m q = 0, 5.5, 40 MeV respectively, and T x = 297(1) MeV, T D = 290(1) MeV for m q = 300 MeV. 
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